Preprint typeset in JHEP style - PAPER VERSION DFTT 27/02 

CERN-TH/2002-191 
IPPP/02/48 
DCPT/02/96 
SHEP-02/24 
March 2003 



Weak Corrections to Three-Jet Production in 

m ■ 

Electron-Positron Annihilations: 1) The Factorisable 
Contributions* 

^ ■ e. Maina 

. Dipartimento di Fisica Teorica - Universita di Torino 

£> \ Via Pietro Giuria 1, 10125 Torino, Italy 

and 

Istituto Nazionale di Fisica Nucleare - Sezione di Torino 
\ Via Pietro Giuria 1, 10125 Torino, Italy 

E-mail: maina@to.infn.it 

(N ■ 

S. Moretti 



in 
o 



C3 



CERN - Theory Division 
CH-1211 Geneva 23, Switzerland 



Qh. and 

' Institute for Particle Physics Phenomenology - University of Durham 

South Road, Durham DEI 3LE, UK 

E-mails: stefano.moretti@cern.ch, stefano.moretti@durham.ac.uk 
j_l ' D.A. Ross 



Department of Physics and Astronomy - University of Southampton 
Highfield, Southampton SO 17 IB J, UK 
Email: dar@hep.phys.soton.ac.uk 

Abstract: We report on the calculation of the factorisable one-loop weak-interaction 
corrections to the initial and final states for three-jet observables in electron-positron an- 
nihilations. We show that such corrections are of a few percent at \fs = Mz- Hence, while 
their impact is not dramatic in the context of LEP1 and SLC, where the total error on the 
measured value of as is larger, at a future Linear Collider, running at the Z mass peak 
(e.g., GigaZ), they ought to be taken into account in the experimental fits, as here the 
uncertainty on the value of the strong coupling constant is expected to be at the 0.1% level 
or even smaller. The calculation has been performed using helicity amplitudes so that it 
can be applied to the case of polarised beams. 
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1. Electroweak corrections at high energies 

Strong (QCD) and Electroweak (EW) interactions are two fundamental forces of Nature, 
the latter in turn unifying weak and electromagnetic (EM) interactions. Together they 
constitute the Standard Model (SM) of particle physics. A clear hierarchy exists between 
the strengths of the two interactions at the energy scales probed by past and present 
high energy particle accelerators (e.g., LEP, SLC, HERA, RHIC and Tevatron) or, indeed, 
at a future generation electron-positron Linear Collider (LC) [1], running with very high 
luminosity at y/s = Mz (the so-called 'GigaZ' stage, where s is the collider CM energy 
squared): QCD forces are stronger than EW ones. This is quantitatively manifest if one 
recalls that the value of the QCD coupling 'constant', as, measured at these machines is 
much larger than the EW one, «EW) typically, by an order of magnitude. 

A peculiar feature distinguishing QCD and EW effects in higher orders is that the 
latter are enhanced by double logarithmic factors, log 2 (-r^-), which, unlike in the former, 



do not cancel for 'infrared-safe' observables [2] - [4]. The origin of these 'double logs' is 
understood. It is due to a lack of cancellation of infrared (both soft and collinear) virtual 
and real emission in higher order contributions. This is in turn a consequence of the 
violation of the Bloch-Nordsieck theorem in non-Abelian theories [5] 1 . The problem is in 
principle present also in QCD. In practice, however, it has no observable consequences, 
because of the final averaging of the colour degrees of freedom of partons, forced by their 
confinement into colourless hadrons. This does not occur in the EW case, where, e.g., the 

'Work supported in part by the U.K. Particle Physics and Astronomy Research Council (PPARC), by the 
European Union (EU) under contract HPRN-CT-2000-00149 and by the Italian Ministero dell'Istruzione, 
dell'Universita e della Ricerca (MIUR) under contract 2001023713.006. 

1 Recently, it has been found that Bloch-Nordsieck violation can also occur in spontaneously broken 
Abelian gauge theories, if the incoming particles are mass eigenstates that do not coincide with gauge 
eigenstates [6]. In the SM this is particularly relevant for incoming longitudinal gauge bosons or Higgs 
scalars [7]. 
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initial state has a non-Abelian charge, dictated by the given collider beam configuration, 
such as in e + e~ collisions. 

These logarithmic corrections are finite (unlike in QCD), as the masses of the EW 
gauge bosons provide a physical cut-off for W and Z emission. Hence, for typical experi- 
mental resolutions, softly and collinearly emitted weak bosons need not be included in the 
production cross-section and one can restrict oneself to the calculation of weak effects orig- 
inating from virtual corrections and affecting a purely hadronic final state. Besides, these 
contributions can be isolated in a gauge-invariant manner from EM effects [4], [8] - [12], 
at least in some simple cases (including the process e + e~ — > j*,Z — ► jets considered here) 
and therefore may or may not be included in the calculation, depending on the observable 
being studied. (See Refs. [8] - [33] for a collection of papers dealing with resummed, one- 
and two-loop EW corrections to various high energy processes.) 

2. One-loop weak effects in three-jets events at leptonic colliders 

It is the aim of our paper to report on the computation of one-loop weak effects entering 
three-jet production in electron-positron annihilation at yfs = Mz via the subprocess 
e + e~ — ► 7*,Z — ► qqg 2 , when the higher order effects arise only from initial or final state 
interactions. These represent the so-called 'factorisable' corrections, i.e., those involving 
loops not connecting the initial leptons to the final quarks, which are the dominant ones at 
yfs = Mz (where the width of the Z resonance provides a natural cut-off for off-shellness 
effects). The remainder, 'non-factorisable' corrections, while being negligible at yfs = Mz, 
are expected to play a quantitatively relevant role as yfs grows larger. (The study of the 
full set of one- loop weak corrections will be the subject of a future publication.) As a whole, 
one-loop weak effects will become comparable to QCD ones at future LCs running at TeV 
energy scales 3 . In contrast, at the Z mass peak, where no logarithmic enhancement occurs, 
one-loop weak effects are expected to appear at the percent level, hence being of limited 
relevance at LEP1 and SLC, where the final error on as is of the same order or larger [36], 
but of crucial importance at a GigaZ stage of a future LC, where the relative accuracy 
of as measurements is expected to be at the 0.1% level or smaller [37]. On the subject 
of higher order QCD effects, it should be mentioned here that a great deal of effort has 
recently been devoted to evaluate two-loop contributions to the three-jet process (albeit, 
only at the amplitude level so far, as there are no numerical results available yet) while 
the one-loop QCD results have been known for some time [38]. Even though a full 0(a|) 
analysis is not yet available, one can reasonably argue that at yfs = Mz the two-loop QCD 
corrections should be comparable to the one-loop weak effects computed here. 

In the case of e + e _ annihilations, the most important QCD quantity to be extracted 
from multi-jet events is precisely as- The confrontation of the measured value of the 

2 See Ref. [34] for the corresponding weak corrections to the Born process e + e~ — > qq and Ref. [35] for 
the ~ m component of those to e + e~ — » qqgg (where m represents the number of light flavours). For 
two- loop results on the former, see [22]. 

3 For example, at one-loop level, in the case of the inclusive cross-section of c + e _ into hadrons, the QCD 
corrections are of 0(2a.), whereas the EW ones are of 0(^f^- log ^-), where s is the collider CM energy 
squared, so that at yfs = 1.5 TeV the former are identical to the latter, of order 9% or so. 
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strong coupling constant with that predicted by the theory through the renormalisation 
group evolution is an important test of the SM or else an indication of new physics, whose 
typical mass scale is larger than the collider energy, but which can manifest itself through 
virtual effects. Jet-shape observables, which offer a handle on non-perturbative QCD effects 
via large power corrections, would be affected as well. 

A further aspect that should be recalled is that weak corrections naturally introduce 
parity-violating effects in jet observables, detectable through asymmetries in the cross- 
section, which are often regarded as an indication of physics beyond the SM. These effects 
are further enhanced if polarisation of the incoming beams is exploited. The option of 
exploiting beam polarisation is one of the strengths of the LC projects. Comparison of 
theoretical predictions involving parity-violation with future experimental data is regarded 
as another powerful tool for confirming or disproving the existence of some beyond the SM 
scenarios, such as those involving right-handed weak currents and/or new massive gauge 
bosons. 

The plan of the rest of the paper is as follows. In the next Section, we describe the 
calculation. Then, in Sect. 4, we present some numerical results. We conclude in Sect. 5. 

3. Calculation 

Since we are considering weak corrections that can be identified via their induced parity- 
violating effects and since we wish to apply our results to the case of polarised electron 

and/or positron beams, it is convenient to work in terms of helicity matrix elements (MEs). 

(G) 

Thus, we define the helicity amplitudes A\ \ 2 a for a gauge boson of type G (hereafter, 
a virtual photon 7* or a Z-boson) of helicity Ai decaying into a gluon with helicity A2, a 
massless quark with helicity a and a massless antiquark with opposite helicity 4 . Since the 
photon is off-shell, it can have a longitudinal polarisation component, so that the helicity 
Ai always takes three values, ±1, 0, for both the 7* and Z gauge vectors 5 , whereas A2 and 
a can only be equal to ±1. 

The general form of these amplitudes may be written as 

<U = n(P2)T^f^v( Pl ), (3.1) 

where p\ and P2 are the momenta of the outgoing antiquark and quark respectively and 
T stands for a sum of strings of Dirac 7— matrices with coefficients, which, beyond tree 
level, involve integrals over loop momenta. Since the helicity a of the fermions is conserved 
the strings must contain an odd number of 7— matrices. Repeated use of the Chisholm 

4 Note that all interactions considered here preserve the helicity along the fermion line, including those in 
which Goldstone bosons appear inside the loop, since these either occur in pairs or involve a mass insertion 
on the fermion line. 

5 These helicities, wherein ±1(0) are(is) transverse(longitudinal), are defined in a frame in which the 
particle is not at rest, so that a fourth possible polarisation in the direction of its four-momentum is 
irrelevant since its contribution vanishes by virtue of current conservation. 
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identity 6 means that F can always be expressed in the form 



r = Ci 7 -pi + C 2 7-P2 + C 3 7-P3 + C 4 VQ 2 ^-n, (3.2) 

where p 3 is the momentum of the outgoing gluon, Q 2 = (p\ + P2 + Pz) 2 is the square 
momentum of the gauge boson, and n is a unit vector normal to the momenta of the jets, 
more precisely: 

V2 pi ■ p 2 pi ■ P3 P2 ■ P3 

The coefficient functions Cj depend on the helicities Ai, A2, a as well as the energy fractions 
x\ and %2 of the antiquark and quark in the final state, i.e., 

OR- 

Xi = ^ (2 = 1,2), (3.4) 
V s 

and on all the couplings and masses of particles that enter into the relevant perturbative 
contribution to the amplitude. 

For massless fermions the MEs of the first two terms of eq. (3.2) vanish, and we are 
left with 

(r\ (1 + 0-7 5 ) (1 + <77 5 ) 

Am M ,o = C 3 u(P2h-P3 { 2 ! v{pi) + C^u(p 2 h-n { 2 } - v(pi), 

= C 3 Q 2 v / (l-x 1 )(l-x 2 ) - iaC 4 Q 2 ^ Xl +X2-l. (3.5) 

The relevant coefficient functions C3 and C4 are scalar quantities and can be projected 
on a graph- by- graph basis using the projections 

Cs = Tv(r 7 .„li^l), (3.6) 

where v is the vector 

(1 - X 2 )pi + (1 - Xi)p 2 ~ (Xi +X 2 - 1)P3 

" : 2Q2(i_ Xl )(i_ X2 ) 

and 

(G) 

At tree level the helicity amplitudes are only functions of xi, X2, the EW couplings g) 
of the (anti)quark of type j (proportional to #w = V^raEW) with «ew = ckem/ sin 2 9w, 
and carrying information on both helicity and flavour of the latter) to the relevant gauge 



6 This identity is only valid in four dimensions. In our case, where we do not have infrared (i.e., soft 
and collinear) divergences, it is a simple matter to isolate the ultraviolet divergent contributions, which are 
proportional to the tree-level MEs, and handle them separately. However, in d dimensions one needs to 
account for the fact that there are 2 d ^ A helicity states for the fermions and (d — 2) for the gauge bosons. 
The method described here will not correctly trap terms proportional to (d — 4) in coefficients of divergent 
integrals. It is probably for this reason that the formalism of Ref. [67] is considerably more cumbersome 
than that presented here. 
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boson and the QCD coupling g$ = s/Airas- Specifically, in case of massless (anti)quarks 
(i.e., m q = 0), we have (here, r a represents a colour matrix): 



with all others being zero. These zero values do not, in general, remain zero in the presence 
of weak corrections and this can lead to a relative enhancement of the latter, in comparison 
to QCD effects at the same order. 

At one-loop level such helicity amplitudes acquire higher order corrections from the 
self-energy insertions on the fermions and gauge bosons shown in Fig. 1, from the vertex 
corrections shown in Fig. 2 7 and from the box diagrams shown in Fig. 3. As we have 
neglected here the masses of the final-state quarks, such higher order corrections depend on 
the ratio Q 2 /My^, where Q 2 is the square momentum of the gauge boson, as well as the EM 
coupling constant oem arid the weak mixing angle s\y = sm6\y (with chew = ckem/s^/)- 
Furthermore, in the case where the final state fermions are 6-quarks, the loops involving 
the exchange of a VF-boson lead to effects of virtual t-quarks, so that the corrections also 
depend on the ratio m^/M/L. (It is only in this case that the graphs involving the exchange 
of the Goldstone bosons associated with the W-boson graphs are relevant.) 

The self-energy and vertex correction graphs contain ultraviolet divergences. These 
have been subtracted using the 'modified' Minimal Subtraction (MS) scheme at the scale 
fi = Mz- Thus the couplings are taken to be those relevant for such a subtraction: e.g., 
the EM coupling, OEM) has been taken to be 1/128 at the above subtraction point. Two 
exceptions to this renormalisation scheme have been the following: 

1. the self-energy insertions on external fermion lines, which have been subtracted on 
mass-shell, so that the external fermion fields create or destroy particle states with 
the correct normalisation; 

2. the mass renormalization of the Z-boson propagator, which has also been carried out 
on mass-shell, so that the Z mass does indeed refer to the physical pole-mass. 

All these graphs are infrared and collinear convergent so that they may be expressed 
in terms of Passarino-Veltman [39] functions which are then evaluated numerically. The 
expressions for each of these diagrams have been calculated using FORM [40] and checked 
by an independent program based on FeynCalc [41]. For the numerical evaluation of the 
scalar integrals we have relied on FF [42]. A further check on our results has been carried 
out by setting the polarisation vector of the photon proportional to its momentum and 

7 Note that we also include self-energy and vertex corrections to the incoming e + e~ — > J* , Z current, 
though we do not show the corresponding graphs. 




(3.8) 



- 5 - 



verifying that in that case the sum of all one- loop diagrams vanishes, as required by gauge 
invariance. The full expressions for the contributions from these graphs are too lengthy to 
be reproduced here. 

In terms of the helicity MEs we define the following "spin-matrix" tensors, only de- 
pending on the polarisation state of the off-shell gauge boson, 



r (GG') 
XX' 



V^4 (G) (a (G ' ] V 

/ -/ A,A? ,cr \- M -X',X 2 ,cr) ' 



(3.9) 



where the (anti)quark and gluon helicities have been summed over. These tensor elements 
are real at tree level, but in general acquire an imaginary part at one loop arising from 
the cuts of the loop integrations above the threshold for the production of the internal 
particles. 

Finally, we define the customary nine form- factors, F±, ...Fg, describing the differential 
structure of a three-jet final state in terms of the above spin-matrix tensors, as follows: 

2 

J 



«EM 
5127T 3 



AA 



Q 2 



(l-Ae 



4s 2 



fi 



L{R) L(R) / f AZ , f ZA 



+ 



where 



1 — A f 



4^: 



Va n z ■ {ir J 
Q 2 (v z 



4sw cw 



(Q2 _ M| + iT z M z ) 



fi 
f. 



GG' 



{ (Q 2 -M 2 z + iT z M z y 

( t-GG' , t-GG' , t-GG' \ 



(i = l,...9)(3.10) 



GG 1 
2 



f GG' 

f GG' 
/4 



n-GG' 
7 0,0 > 

-2Ke (Tff 



n-GG' 
2 -l,-l 



GG' 



GG' 



f 
f GG 

h 



GG' 



GG' 



-^(r^' + r^'), 

-2^9^(7^+7^), 



(3.11) 



with 7/^^ the weak correction factor to the coupling of the left(right)-handed electron to 
the gauge boson G and A e the helicity of the incoming electron beam (assumed always to 
be of opposite helicity to the incoming positron beam). 

Up to an overall constant, these form- factors are the same as those introduced, e.g., 
in Refs. [43,44]. The last three (F7, ...Fg) can arise for the first time at the one-loop level, 
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since they are proportional to the imaginary parts of the spin-matrix. Besides, F3, Fq 
and F7 vanish in the parity-conserving limit and can therefore be used as probes of weak 
interaction contributions to three-jet production. (Moreover, F3 and Fq would be exactly 
zero at tree level if the leading order process were only mediated by virtual photons.) 

These form-factors further generate the double differential cross-section for three-jet 
production in terms of some event shape variable, S, which is in turn related to x±, X2 by 
some function, s, i.e., S = s(xi,X2), and of the polar and azimuthal angles, a, (3, between, 
e.g., the incoming electron beam and the antiquark jet, by 8 

J dxidx 2 S(S - s(x 1 ,x 2 )) [(2 - sin 2 a) F x + (1 - 3 cos 2 a) F 2 

+ A e cos a F3 + sin 2a cos (3 F4 + sin 2 a cos 2f3 F5 + A e sin a cos (3 Fq 
+ sin 2a sin (3 F-j + sin 2 a sin 2/3 Fg + A e sin a sin (3 Fg] . (3-12) 

Note that upon integrating over the antiquark jet angle relative to the electron beam, only 
the form-factor iq survives. 

In general, it is not possible to distinguish between quark, antiquark and gluon jets, 
although the above expression can easily be adapted such that the angles a, (3 refer to 
the leading jet. However, (anti)quark jets can be recognised when they originate from 
primary 6-(anti)quarks, thanks to rather efficient flavour tagging techniques (such as /j,- 
vertex devices) . We will therefore consider the numerical results for such a case separately. 

4. Numerical results 

The processes considered here are the following: 

e + e~ — ► 7*, Z — ► qqg (all flavours), (4-1) 

when no assumption is made on the flavour content of the final state, so that a summation 
will be performed over q = d,u, c, s, 6-quarks, and also 

e+e- ^ 7 *,Z^bbg, (4.2) 

limited to the case of bottom quarks only in the final state. As already intimated, all 
quarks in the final state of (4.1)-(4.2) are taken as massless 9 . In contrast, the top quark 
entering the loops in both reactions has been assumed to have the mass mt = 175 GeV. 
The Z mass used was Mz = 91.19 GeV and was related to the VF-mass, M\y, via the SM 
formula Mw = Mz cosOw, where sin 2 9w = 0.232. (Corresponding widths were Tz = 2.5 
GeV and IV = 2.08 GeV.) For as we have used the two-loop expression for A^i -4 -* = 200 
MeV throughout (yielding, a s (M z ) = 0.11). 

8 A qualitative difference between the expressions of the form-factors, Fi (i = 1, ...9), used here and 
those of Refs. [43,44] is that we do not include the sign of the axial vector coupling of the electron to the 
exchanged gauge boson in our definitions. In this way the difference between the differential cross-sections 
for left- and right-handed polarised electron beams is manifest in eq. (3.12). 

9 Mass effects in e + e~ — » 7*, Z — » bbg have been studied in [45] and [46]. 
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We systematically neglect higher order effects from EM radiation, including those due 
to Initial State Radiation (ISR) or beamstrahlung. In fact, although these are known to be 
non- negligible (especially at LC energies), we expect them to have a similar effect on both 
the tree-level and one-loop descriptions, hence being irrelevant for our purpose. In this 
context, we should like to elaborate further on the purely EM corrections to the final state 
of processes (4.1)— (4.2). Those to the form-factor F\ have already been calculated, since 
they can be extracted from the Abelian part of the NLO-QCD corrections (see [38]) by 
replacing Cf by unity and as by oem- As was pointed out in Ref. [47], these corrections are 
dominated by a term ~ multiplying the tree- level cross-section. This contribution 

is ~ 1% and is independent of the jet event shape. A further correction, associated with 
the Sudakov form-factor, acts in the negative direction and is subdominant away from the 
two-jet region (i.e., up to values of ~ 0.95 for the Thrust, see below for its definition). There 
is no reason to believe that these EM corrections would be enhanced for other form-factors. 

It is common in the specialised literature to define the n-jet fraction R n {y) as 

Rn(v) = (4.3) 

where y is a suitable variable quantifying the space-time separation among hadronic objects 
and with Co identifying the (energy-dependent) Born cross-section for e + e~ — > qq. 

For the choice [i = sj~s of the renormalisation scale, one can conveniently write the 
three-jet fraction in the following form: 

MV) = (g) A (v) + (5) 2 B(y) + (4.4) 

where the coupling constant as and the functions A(y) and B{y) are defined in the MS 
scheme. An experimental fit of the R n {y) jet fractions to the corresponding theoretical 
prediction is a powerful way of determining as from multi-jet rates. 

Through order 0(as) processes (4.1)-(4.2) are the leading order (LO) perturbative 
contributions to the corresponding three-jet cross-section 10 , as defined via eqs. (4.3)-(4.4). 
The LO terms, however, receive higher order corrections from both QCD and EW interac- 
tions and we are concerned here with the next-to- leading order (NLO) ones only. Whereas 
at LO all the contributions to the three-jet cross-section come from the tree-level parton 
process e + e~ — > qqg (which contributes to the A{y) function above), at NLO the QCD con- 
tributions to the three-jet rate (hereafter, denoted by NLO-QCD) are due to two sources. 
First, the real emission diagrams for the processes e + e~ — ► qqgg and e + e~ — > qqQQ, 
in which one of the partons is 'unresolved'. This can happen when one has either two 
collinear partons within one jet or one soft parton outside the jet. Both these contri- 
butions are (in general, positively) divergent. Thanks to the Bloch-Nordsieck [48] and 
Kinoshita-Lee-Nauenberg [49] theorems, these collinear and soft singularities are cancelled 
at the same order in as by the divergent contributions (generally negative) provided by the 
second source, namely, the virtual loop graphs. Therefore, after renormalising the coupling 
constant as, a finite three-jet cross-section is obtained and the function B(y) accounts for 



'Hereafter, perturbative contributions are refereed to relatively to the 0(a| M ) two-jet rate. 
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the above-mentioned three- and four-parton QCD contributions 11 . While the EM compo- 
nent of the EW corrections may be treated on the same footing as the QCD one (with 
the additional photon playing the role of a second gluon) , the weak corrections of interest 
(hereafter, labelled as NLO-W) only contribute to three-parton final states. Hence, in order 
to account for the latter, it will suffice to make the replacement 

A(y) - A(y)+A w (y) (4.5) 

in eq. (4.4). 

The decision as to whether two hadronic objects are unresolved or otherwise is usually 
taken through the application onto the hadronic final state of a so-called 'jet clustering 
algorithm', wherein the number of clusters 12 is reduced one at a time by combining the 
two most (in some sense) nearby ones. The joining procedure is stopped by testing against 
some criterion and the final clusters are called jets. 

As jet clustering schemes 13 , we have used a selection of the 'binary' ones, in which 
only two objects are clustered together at any step. Given two clusters labelled as i and 
j, the measure of their 'distance' is normally denoted by yij and the minimal separation 
allowed by y cut . The algorithms are the following: the JADE (J) one [51], which uses as a 
measure of separation the quantity 



s 



the Durham (D) [52] and the Cambridge (C) [53] ones, both using 



11 



D c 2mm(Ei,Ej)(l-cos0 ij ) 

'hi = Vij = - B ; ( 4 -7) 



the Geneva (G) one [55], for which one has 

G _ 8 £;£,-(!- cos %) 



1tJ 9 (Ei + EjY 



(4.8) 



In eqs. (4.6)— (4.8), E{ and Ej are the energies and 0y the angular separation of any 
pair ij of clusters in the final state. The choice of these particular schemes has a simple 
motivation. The D and C ones are different versions of 'transverse-momentum' based 
algorithms, whereas the J and G ones use an 'invariant-mass' measure (see [54] for a 
review). In fact, these two categories are those that have so far been employed most 
in phenomenological studies of jet physics in electron-positron collisions, with the former 



11 In order to calculate these, we make use here of a program based on Ref. [50]. 

12 Here and in the following, the word 'cluster' refers to hadrons or calorimeter cells in the real experi- 
mental case, to partons in the theoretical perturbative calculations and also to intermediate jets during the 
clustering procedure. 

13 We acknowledge here the well admitted abuse of referring to the various jet 'finders' both as algorithms 
and as schemes, since the last term was originally intended to identify the composition law of four-momenta 
when pairing two clusters: in our case, = pf + p 1 ^ . 

14 The Cambridge algorithm in fact only modifies the clustering procedure of the Durham jet finder and 
the two implementations coincide for n < 3 parton final states. 
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gradually overshadowing the latter, thanks to their reduced scale dependence in higher 
order QCD (e.g., in the case of the 0(a§) three- [54] - [56] and 0(a\) four-jet rates [57]) 
and to smaller hadronisation effects in the same contexts (see Refs. [54,55]). 

Fig. 4 displays the A(y), —Aw(y) and B(y) coefficients entering eqs. (4.4)-(4.5), as 
a function of y(= y cut ) for the four above jet algorithms at yfs = Mz 15 - A comparison 
between A(y) and Ayj(y) reveals that the NLO-W corrections are negative and remain 
indeed at the percent level, i.e., of order 2 °™ without any logarithmic enhancement (since 
^/s ~ M\y, Mz)- They give rise to corrections to 03(2/) of -1%, and thus are generally 
much smaller than the NLO-QCD ones. In this context, no systematic difference is seen 
with respect to the choice of jet clustering algorithm, over the typical range of application 
of the latter at yfs = Mz (say, y cut ^ 0.005 for D, C and y cut ^ 0.01 for G, J). 

As already mentioned, it should now be recalled that jets originating from 6-quarks can 
efficiently be distinguished from light-quark jets. Besides, the 6-quark component of the full 
three-jet sample is the only one sensitive to i-quark loops in all diagrams of Figs. 1-3, hence 
one may expect somewhat different effects from weak corrections to process (4.2) than to 
(4.1) (the residual dependence on the Zqq couplings is also different). This is confirmed 
by Fig. 5, where we present the total cross section at y/s = Mz for e + e~ — ► 7*,^ — ► bbg 
as obtained at LO and NLO-W, for our usual choice of jet clustering algorithms and 
separations. A close inspection of the plots reveals that NLO-W effects can reach the 
~ —2.0% level or so. 

In view of these percent effects being well above the error estimate expected at a future 
high-luminosity LC running at the Z pole, it is then worthwhile to further consider the 
effects of NLO-W corrections to some other 'infrared-safe' jet observables typically used 
in the determination of as, the so-called 'shape variables' [58]. A representative quantity 
in this respect is the Thrust (T) distribution [59]. This is defined as the sum of the 
longitudinal momenta relative to the (Thrust) axis rax chosen to maximise this sum, i.e.: 

T = max^^pH, (4.9) 

Ei \Pi\ 

where i runs over all final state clusters. This quantity is identically one at Born level, 
getting the first non-trivial contribution through 0(«s) from events of the type (4.1)— (4.2). 
Also notice that any other higher order contribution will affect this observable. Through 
O(ag), for the choice (i = yfs of the renormalisation scale, the T distribution can be 
parametrised in the following form: 

<'-^;Hs) AT(T) + (s) V(T) - 

Again, the replacement 

A T (T)^A T {T) + A^(T) (4.11) 
accounts for the inclusion of the NLO-W contributions. 



15 Notice that A(y) and Aw(y) for the C scheme are identical to those for the D one (recall the previous 
footnote). 
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We plot the terms (ff ) A T (T), (§a) A^(T) and (g) 2 5 T (T) in Fig. 6, always at 
= Mz", alongside the relative rates of the NLO-QCD and NLO-W terms with respect 
to the LO contribution. Here, it can be seen that the NLO-W effects can reach the level of 
— 1% or so and that they are fairly constant for 0.7 ^ T ^ 1. For the case of 6-quarks only, 
similarly to what seen already for the inclusive rates, the NLO-W corrections are larger, 
as they can reach the —1.6% level. 

The ability to polarise electron (and possibly, positron) beams joined with the high lu- 
minosity available render future LCs a privileged environment in which to test the structure 
of hadronic samples. As noted earlier, differential spectra may well carry the distinctive 
hallmark of some new and heavy strongly interactive particles (such as squarks and gluinos 
in Supersymmetry) , whose rest mass is too large for these to be produced in pairs as real 
states but that may enter as virtual objects into multi-jet events. Similar effects may how- 
ever also be induced by the NLO-W corrections tackled here. Both could well be isolated 
in one or more of the nine form- factors given in eq. (3.10). As already intimated in the 
previous Section, F-j to F$ are identically zero at LO 16 , even prior to any integration in a, (3 
and/or averaging over the e + e~ helicities. Besides, Fj would remain zero unless correc- 
tions involve parity- violating interactions whereas F% , Fq and Fg would not contribute for 
unpolarised beams. As for F±, ...Fq, we should mention that the NLO-W corrections to the 
corresponding tree-level distributions were found to be < 1.2(2.0)% for left-(right-)handed 
incoming electrons 17 . 

Observables where such effects would immediately be evident are what we call the 
'unintegrated' (or 'oriented') Thrust distributions associated to each of the form-factors 
in eq. (3.12) (wherein S = T). In Fig. 7, we present the Fi = i^(T) terms appearing in 
that expression, each divided by ctq (for consistency with the previous plots), alongside 
the absolute value of the relative size of the NLO-W corrections with respect to the LO 
case, for the form-factors Fi,...Fq, which are non-zero at the Born level. For the latter, 
NLO-W corrections can be either positive or negative, depending on the form-factor being 
considered, and can be as large as ±4% or so (in the case of F3 and F$). 

5. Conclusions 

On the basis of our numerical findings in the previous Section, we should like to conclude 
as follows. 

• At ^/s = Mz, the size of the NLO-W corrections to three-jet rates is rather small, of 
order percent or so, hence confirming that determinations of as at LEP1 and SLC are 
stable in this respect and that the SM background to parity-violating effects possibly 
induced by new physics is well under control. In contrast, NLO-W effects ought to 
be included in the case of future high-luminosity LCs running at the Z pole, such as 

16 This is strictly true only for massless quarks, as, for m q =fc 0, Ref. [43] has shown that Fg becomes 
non-zero. 

17 For reasons of space we refrain from presenting here the NLO-W dependence of Fi, ...Fc, in term of x\ 
and X2- The files can be requested from the authors. 
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GigaZ, where the accuracy of as measurements from jet rates is expect to reach the 
0.1% level. 

• Exclusive observables in three-jet events are also affected by similar NLO-W effects: 
e.g., the Thrust distribution, as representative of the so-called 'infrared safe' quanti- 
ties. The experimental error expected at LCs in the determination of as from such 
quantities at y/s = Mz is again of the order of 0.1% (or even smaller), so that the 
inclusion of NLO-W effects in the corresponding theory predictions is then manda- 
tory. 

• Effects from NLO-W corrections are somewhat larger in the case of 6-quarks in the 
final state, in comparison to the case in which all flavours are included in the hadronic 
sample, because of the presence of the top quark in the one-loop virtual contributions. 

• Since the exploitation of beam polarisation effects will be a key feature of experimen- 
tal analyses of hadronic events at future LCs, we have computed the full differential 
structure of three-jet processes in the presence of polarised electrons and positrons, 
in terms of the energy fractions of the two leading jets and of two angles describing 
the final state orientation. The cross-sections were then parametrised by means of 
nine independent form-factors, the latter presented as a function of Thrust at fixed 
angles. Three of these form-factors carry parity- violating effects which cannot then 
receive contributions from ordinary QCD. For the two that are non-zero at LO, the 
NLO-W corrections were found as large as 4%. Such higher-order weak effects should 
appropriately be subtracted from hadronic samples in the search for physics beyond 
the SM. 

• All our results were presented for the case of the factorisable NLO-W effects, i.e., 
for corrections to the initial and final states only. Whereas these should be sufficient 
to describe adequately the phenomenology of three-jet events at LEP1, SLC and 
GigaZ energies, at TeV energy scales one expects comparable effects due to the non- 
factorisable corrections, in which weak gauge bosons connect via one- loop diagrams 
electrons and positrons to quarks and antiquarks. Their computation is currently in 
progress and we will report on it in due course. 
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Figure 1: Self-energy insertion graphs. The shaded blob on the incoming wavy line represents all 
the contributions to the gauge boson self-energy and is dependent on the Higgs mass (hereafter, we 
will use Mh = 115 GeV for the latter). In this and all subsequent figures the graphs in which the 
exchanged gauge boson is a V^-boson is accompanied by corresponding graphs in which the W-boson 
is replaced by its corresponding Goldstone boson. Since the Yukawa couplings are proportional to 
the fermion masses, such graphs are only significant in the case of 6-quark jets. There is a similar 
set of diagrams in which the direction of the fermion line is reversed. 
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Figure 2: Vertex correction graphs. Again, same considerations as in the previous figure apply 
for the case of Goldstone bosons and there is a similar set of graphs in which the direction of the 
fermion line is reversed 




Figure 3: Box graphs. Again, same considerations as in the previous two figures apply for the case 
of Goldstone bosons. Here, the first graph is accompanied by a similar graph with the direction 
of the fermion line reversed whereas for the second graph this reversal does not lead to a distinct 
Feynman diagram. 
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Figure 4: The A(y), —Aw and B(y) coefficient functions of eqs. (4.4)-(4.5) for the Cambridge, 
Durham, Geneva and Jade jet clustering algorithms, at y/s — Mz- (Notice that the ~ Aw term 
has been plotted with opposite sign for better presentation.) 
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Figure 5: The total cross section for process (4.2) at LO and NLO-W for the Cambridge, Durham, 
Geneva and Jade jet clustering algorithms, at y/s — Mz- (Notice that the NLO-W results have 
been plotted with opposite sign for better presentation.) 
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Figure 6: The LO, NLO-QCD and NLO-W contributions to the coefficient functions entering the 
integrated Thrust distribution, see eq. (4.10), for process (4.1) (top) and the relative size of the 
two NLO corrections (bottom), at y/s = Mz- The correction for the case of &-quarks only is also 
presented, relative to the LO results for process (4.2). (Notice that the ~ Aw terms have been 
plotted with opposite sign and multiplied by hundred for better presentation.) 
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e + e~ -> 3 jets at NLO-W 
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Figure 7: The unintegrated Thrust distributions for the nine component of the cross-section 
associated to the form-factors in eq. (3.10) for the NLO-W process (4.1) (top and middle) and the 
relative size of the six components which are non-zero at LO (bottom), at i/s = Mz- Labels are 
as follows: (top) F 4 (solid), F 3 (short-dashed); (middle) F 9 (solid), F 8 (short-dashed), F7 (dotted), 
F 6 (dot-dashed), F 5 (dashed), F 2 (fine-dotted), F\ (long-dashed); (bottom) F e (solid), F 5 (short- 
dashed), F 4 (dotted), F 3 (dot-dashed), F 2 (dashed), F\ (fine-dotted). 
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